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Abstract 

We find the best asymptotic lower bounds for the coefficient of the leading term of 
the L\ norm of the two-dimensional (axis-parallel) discrepancy that can be obtained 
by K. Roth's orthogonal function method among a large class of test functions. We 
use methods of combinatorics, probability, complex and harmonic analysis. 

1 Introduction 

Let P be a finite set in the two-dimensional unit square (P C [0, l] 2 ). Define the 
two-dimensional (axis-parallel) discrepancy function as follows: 

D(x, y) = D P (x, y) = fj(P n ([0, x] x [0, y})) - (t(P) • xy, (x, y) £ [0, l] 2 , 

where fl(Q) denotes the number of elements (cardinality) of a set Q. The discrepancy 
function measures the difference between the actual number of points of the set P in 
an axis-parallel rectangle [0, x] x [0, y] and the "expected" number of points in that 
rectangle. Thus, the discrepancy function quantifies the "closeness" of the distribution 
generated by a finite set P to the uniform distribution. 

The Lqo norm of the discrepancy function (also known as the star-discrepancy) - 
I|PpIIloc([o,i] 2 ) has been studied for a while and is known to be related to numerical 
integration, metric entropy, the small ball inequality etc. We refer the reader to |4J,[6J 
or [9] for a comprehensive discussion of the star-discrepancy. 

In this paper we will be dealing with lower bounds for the L\ norm of the discrepancy 
function. Quoting a recent essay on the discrepancy theory ([1], 2013): "The other 
endpoint of the L p scale, p = 1, is not any less (and perhaps even more) difficult than the 
star-discrepancy estimates. The only information that is available is the two-dimensional 
inequality (proved in the paper of Halasz [3j,1981))". 

When discussing the L\ norm of the discrepancy function, it is convenient to define: 
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djy = — ? • inf / / \Dp(x, y)\ dxdy. (1-1) 

V^n m=Nj Jo 

It can be deduced from H. Davenport's argument ([2], 1956) that: 

limsupdjv < +oo. (1-2) 
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On the other hand (see remark |4T3"|) G. Halasz ([3], 1981) proved the following inequality: 

liminf d N > c H > 0. (1.3) 

N— >oo 

A careful analysis of his proof shows that: 

c H = 7= « 0.00039. 

1152(v / e + l)v / hi2 



The inequalities (|1.2p and (|1.3p show that d/v is the coefficient of the leading term of the 
L\ norm of the two-dimensional (axis-parallel) discrepancy function. 

In this paper we improve the asymptotic lower bounds for d^, namely: 

Theorem 1.4. For the coefficient of the leading term of the L\ norm of the two- 
dimensional discrepancy defined by (jl.ip the following asymptotic lower bounds hold: 

liminf d N > -. ra 0.00854 

n->°c 256Veln2 

and 

limsupdjv > : ~ 0.01138. 

Af^oo 64\/eln2 

For comparison, the best-known lower estimates of the L2 norm of the discrepancy 
were obtained in [8]. 

Remark 1.5. In jllOf it was proved that: 

liminf -== ■ inf ||L> P || L , ao 1,2-, < 0.17905. 

Numerical experiments for the sets introduced in fTB/ make the following estimate plau- 
sible: 

lim inf , inf ll-Dpll Tn( i n 1121 < 0.17601. 



Remark 1.6. Comparison of remark \1.5\ with theorem \1.4\ shows that for sets with 
minimal L2 discrepancy their L2 discrepancy is close to the L\ discrepancy. 

In this paper we also show that theorem 11.41 provides the best bounds that can be 
obtained by K. Roth's orthogonal function method for a large class of test functions (see 
theorems 14.71 and IT. ip . 

Theorem 11.41 as well as the fact that the best known lower bounds for the lim sup 
and lim inf of the coefficient of the leading term of the L2 norm of the two-dimensional 
discrepancy are also different (this fact can be derived from [8]), lead us to formulate 
the following conjecture: 

Conjecture 1.7. The following strict inequality is true: 

lim sup dj^ > lim inf djy 
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2 Auxiliary Functions 

The existing proofs of lower bounds for discrepancy are based on certain construction 
of functions called Roth's auxiliary functions (see [1]). In this section we will modify 
that construction so that now the auxiliary functions do not take the value zero, but 
from all other perspectives they are as good as Roth's original auxiliary functions. The 
fact that they do not take the value zero will play a crucial role later on. 

For a finite set P C [0, l] 2 denote by N = tf(P) the number of elements of the set P. 
Let n be the only integer satisfying the condition: 

< 2(t(P) = 2N < 2 n . (2.1) 

For every i = 0, 1, . . . , n we will define E{ to be a family of dyadic rectangles in [0, l] 2 : 

+oo 

Ei = |J El 

1=0 

where the sets E\ will be defined recursively in the next paragraph, and we will define 
Ci to be a family of dyadic rectangles in [0, l] 2 : 

Ci = [J c\, 

1=0 

where the sets C\ will be defined recursively in the next paragraph. 
Recursive definition. 

Step 0. In this step we define the families E® and Cf. Consider the dyadic rectangles 
R = R x x R y € [0, l] 2 whose sidelengths satisfy \R X \ = 2~\ \R y \ = 2 i ~ n . Denote by Cf 
the family of those sets R that contain at least a point of the set P and by Ef the family 
of those sets R that do not contain any point of the set P. We have: 

tl(Cf) + t)(P l ) = 2" 

and 

ft(CP) < N. 

Step 1. Partition each of the sets R G Cf into 2 2 ( n+1 ) dyadic rectangles which have 
the form P 1 = P* x R l y , where \R l x \ = \R x \/2 n+1 , |H£| = \R x \/2 n+1 . Define Ej to be the 
family of those sets P 1 which do not contain any point of the set P and define C\ to be 
the family of those sets P 1 which contain at least a point of the set P. We have: 

i{C}) < N. 
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Step l + l. In step I + 1 partition each of the sets R l = R l x x R l y G C\ into 2 2 ( n+1 ) 
dyadic rectangles which have the form R l+1 = R l x +1 x where 

| = \R l x \/2 n+1 , \R y +1 \ = \R l y \/2 n+1 . 

Define E^ +1 to be the family of those sets R l+1 which do not contain any point of the 
set P and define C 4 - +1 to be the family of those sets R 1 which contain at least a point of 
the set P. We have: 

< N. 

Definition 2.2. Let R = R x x R y be a dyadic rectangle in [0, l] 2 . By Hr we will denote 
the -Lqo normalized Haar function adapted to the dyadic rectangle R, i.e. 

hn(x,y) = h Rx (x) ■ h Ry (y), 

where the function h Rx equals 1 on the left half of the dyadic interval R x , —1 on the 
right half of the dyadic interval R x and outside of the dyadic interval R x . Similarly, 
the function h Ry equals 1 on the left half of the dyadic interval R y , —1 on the right half 
of the dyadic interval R y and outside of the dyadic interval R y . 

Definition 2.3. Define the auxiliary functions fi by: 

ReEi 

where Iir is the normalized Haar function adapted to the rectangle R (see definition 
Remark 2.5. The following function: 

ReE° 

is the original Roth 's auxiliary function. 
Lemma 2.6. The set 

M=[0,1] 2 \ |J R 

has measure zero. 
Proof. 

\R\ = tt(C^)2- n - 2 ( n+1 )' = at 2 -"- 2 ("+ 1 )' _). 0, as / -»• +oo. 

ReCj 

□ 

Remark 2.7. Due to lemma [275[ the function fi is defined almost everywhere on [0, l] 2 . 
For definiteness, we will additionally take fi to be equal to 1 on the remaining set of 
measure zero. 
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We now investigate the properties of the functions /j. 
Lemma 2.8. 

h , [ ,l] 2 -^{-l, 1} 

or equivalently 

\fi(x,y)\ = l. 

Proof. The proof follows from the definition of the functions /, (see definition 12, 3p , from 
the remark \2. 71 and the fact that the elements of Ei are disjoint. □ 

Lemma 2.9. 

»1 r l N . 2 -2n 

Dpf t <-{T-N) . (2.10) 

o Jo 16 

Proof. We use the following well-known fact about the interrelation between discrepancy 
and Haar functions (see e.g. |6j): if the coordinate rectangle R does not contain any 
point of the finite set P then: 



Jo 
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Let Ef be as in the construction of the function /j, then we have: 

1 "»w,<e [ [D P . h ^-j:^f<-^-N)^§ 2 " 



JO r^^^nJO JO 



Lemma 2.11. Zei < ii < Z2 < • • • < i p < n. T/ien 

»1 /•! 



□ 



and 

Proof. Let 
Denote: 



/ / fiifi2 ■ ■ ■ fip — 
JO JO 

/ / D P ■ fhfi 2 ■ ■ ■ fip 

Jo Jo 



»1 /-l 

6) 



■ ■ N2~ n 
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-Rii — Phi,x X -^12 — ^12, x -^12 ,3/? ■ • • > -^ip — Ri p ,x X Ri p , 

For definiteness assume: 

ip I _ Q-ii-(n+l)ii ip | _ 911— n— (n+l)Zi 

I -'Ml, a: I — * > l-'Mijj/l — ^ j 

|p | _ o— 12 — (n+l)'2 |p | _ o*2— n— (n+l)J2 

| i >'i2,2:| — * ) |- fL l2 ) 2/l — J 
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IP I 9— ip — (n+l)Z p |p 1 ni p —n—(n+l)lp 

\ I Mp,x\ — * 1 \ ri 'ip,y\ — * 

Note that the repeated subdivisions into 2 2 ( n+1 ) dyadic rectangles in the construction in 
section [2] guarantees that the lengths of the dyadic intervals Ri 1>x , ■ ■ ■ , Ri p ,x are mutually 
distinct, similarly the lengths of the dyadic intervals Ri lt y, ■ ■ ■ , Ri p ,y are mutually distinct. 
Hence the product: 

hn x h R2 . . . hn p 

is either identically zero or (up to a sign) the Haar function on the intersection: 



r = n Ri. 



t=i 



Thus the first claim of the theorem follows. 

We now estimate the sidelengths of the dyadic rectangle R. Write R = R x x R y then: 



and therefore 
Now denote: 



E 



\R X \ < < 2 n ~ n , 

\R\ = \R X x R y \ < 2- l p +il - n . 
p 

R = P| Rit '■ Rii e Eii > ■ ■ ■ ) Rip G Eip 

t=l 



(2.12) 



Note that the elements of E are mutually disjoint (as the elements of each Ei t are 
mutually disjoint). We also have: 



fhfi 2 ■■■U v — 



RaE 



Using (|2TT2l we get: 



/ / D p .f h f i2 ...f ip <J2 [ [ Dp 
Jo Jo ReE Jo Jo 



■hi 



v-^ iV|i?| 2 N2- l p +ll -' n v-^ . „. , , N2 
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ReE 
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□ 
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3 Combinatorial properties of the auxiliary functions 

We now explore the combinatorial properties of the auxiliary functions /, defined bv !2.3l 
We adapt the following notation from the theory of symmetric functions (see e.g. [5]): 

n 
i=0 

e 2 = &2 = fiifi2i 
0<ii<i2<n 

63 = 63 = fh fi2 fi-3 > 

0<ii<i2<«3<n 

and so on. 

Let the number k be odd (k = 1, 3, 5, . . . ). Then using the fact that ff = 1 (see 
lemma [2T8]) we can write: 

(/o + /! + •■■ + f n ) k = A?(*)e? + + Ag(fc)eR +■■■ (3-2) 

(here we assume that A™{k) = once p > min(/c, n + 1)). 
Remark 3.3. For example we compute: 

(/o + /i) 3 = /o 3 + 3/o/i 2 + 3/ 2 /i + A 3 = 4/o + 4/ l5 

so that A\(3) = 4. 

Remark 3.4. It is important to point out that the coefficient of the term of 

{h + fl + --- + fnf 

that is linear with respect to fo + fi + — - + f n is A™{k). 

4 Test Functions 

We start this section by sketching Roth's orthogonal function method (as modified 
by Halasz |3j). 

Let P C [0, l] 2 be a finite set and let /o, /1, • • • , /n be the corresponding functions 
defined bv !2.3l Let H n = H n (zo, z\, . . . , z n ) be an odd, symmetric function that is entire 
in each of its n + 1 variables. 

By LIN(H n ) we will denote the coefficient of the term of the composite function 
H n (fo, /1, . . . , f n ) which is linear with respect to Yli=o /«• For example, according to the 
remarkESl LIN{(f + ■ ■ ■ + f n ) k ) = A?(k). 

Assume that 

\\H n {z ,Z 1 ,...,Z n )\\ Loo ( [ _ 1A]n +l ) = 1. (4.1) 
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Then we can estimate the L\ norm of the discrepancy function Dp{x, y) from below 
as follows: 



\Dp\\i > 



C f 1 D P -H n (foJi,---Jn) > LIN(H n ) V f 1 C D P f, 
Jo Jo i=0 Jo Jo 



+ Error, 



where in Roth's orthogonal function method it is further proved that the term Error is 
negligible (estimating the term Error from above by means of lemma [2.1ip . Then lemma 
12.91 is applied to obtain the following lower bound for the L\ norm of the discrepancy 
function: 

N2~ 2n 

\\D P \\i >LIN(H n ){n + l)(2 n - N)—— + Error, (4.2) 

lb 

where the integers N and n are defined by (|2.ip . 

Informally, Roth's orthogonal function method (as modified by Halasz |3j) consists 
of finding appropriate functions H n (called test functions) for which LIN(H n ) grows as 
fast as possible, whereas the error term Error is negligible, so that an appropriate lower 
estimate for ||.Dp||i follows. 

Remark 4.3. In his proof G. Halasz f3jj used the following test function: 

where the auxiliary functions /• are as in remark \2. 5\ and 7 > is a certain small 
constant. This test function involves complex numbers. The choice of the test function 
G was motivated by Riesz products. We will use a different approach by choosing our 
test function to be a solution of an extremal problem in Fourier analysis (see theorem 

vn\i. 



Remark 4.4. The rest of the paper is organized as follows. We will first show that the 
problem of finding appropriate test functions H n reduces to finding a certain function T 
(see theorem \4-7\ ), so that the quantity corresponding to LIN{H n ) is: 

LIN n {T) = LIN (t ( fo + "^- +f A \ . (4.5) 



Then we will derive a precise formula for LIN n (T) (see section\B^. After that we will 
find the function T that maximizes the quantity: 

lim sup y/n ■ LIN n (T) (4.6) 



in a certain natural class of functions (see theorem \ 7. l\ j and show that for that function 
the error term Error appearing in (|4.2p can be neglected (see section^). 

Theorem 4.7. Without loss of generality one can assume that the test functions H n 
have the following form: 



H n {fo, /!,•••,/« 



fo + fl + ■ ■ ■ + fn 



n 
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Proof. According to the "Main theorem on symmetric functions" (see e.g. [3]) the sym- 
metric function H n (fo, fi,-..,f n ) can be written as a function depending only on the 
functions (introduced by (|3.ip ). We borrow the following notation from the theory of 
symmetric functions: 

0<i<n 

Recalling Newton's identities: 

k 

ke k = yii-lf^ek-iPi, 
i=l 

we obtain that a symmetric function can be written as an analytic function depending 
only on the functions pi. Taking into account that ff = 1 (see lemma l2?8|) we have: 

Pk = Pi, 1 < k < n + 1, k is odd, 

and 

Pk = n + 1 = Po, 1 < k < n + 1, k is even. 

Thus, the function H n depends only on po anci Pi, where po is a constant. Recalling that 
H n is odd, it turns out that H n has to depend only on: 

Pi = /oH h/ n , 

i.e. we can write the function H n in the following form: 

H n (fo, fl, . . . , fn) = F n (f + • • • + /„). 
Introduce functions T n as follows: 



. n 

Then the condition ()4.ip becomes: 

sup \T n {x)\ = 1. 
(-v^.v 7 ") 

Recall that {T n (z)} c £ = i is a family of odd, entire functions. Either it is normal or it 
goes locally uniformly to 0, or to oo. The latter two cases are not interesting, so let us 
assume that T n has a subsequence that tends to a certain odd, entire function T locally 
uniformly in the complex plane. For clarity of notation we will assume that T n itself 
tends to T. 

Note that the functions fa are non-correlated on the unit square and take the values 
1 and —1 with equal probability, thus: 

/o + ■■■ + /„ . jv (0j i). 



■n 
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So that the conditions: 

sup \T n (x)\ = 1 

imposed on the functions T n induce the following condition on the limiting function T : 

sup \T(x)\ = 1, 

— oo<x<+oo 

and instead of discussing the functions: 

7o + ■■■ + /« 



Tn 
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we can discuss the functions: 



rp i fo H \- fn 



n 

□ 



5 Generating function of A^ik) 



In this section we derive the formula (|5.ip for the generating function of the sequence 
{A r l(k)} k= i i 3 ! .„ introduced by the identity (f3T2|) . 

Denote by B n the binomial distribution, i.e. the distribution of the sum of n inde- 
pendent random variables which take values 1 and with probability 1/2. Denote by 
S n the distribution of n independent random variables which take values 1 and — 1 with 
probability 1/2. Then 

S n = 1B n n. 

Denote by Mb„ the moment-generating function of the binomial distribution B n , that 
is: 

oo 1 

M Bn {t) = E(e tB ") =J2u E 

where E{X) is the mathematical expectation of a random variable with distribution X. 
We know that: 

i , KtY 



As a consequence we get: 

M 2 B n -n+l(t) = E f e (2^-n+l)A = e (-' n + 1 ) t M Bn (2t) 



Let xq, x±, . . . , x n be independent random variables, where xq = 1 and xi,...,x r 
take values 1 and —1 with probability 1/2. Let k be an odd number. We have: 

E ((x + ■ ■ ■ + x n ) k ) =A n l {k). 
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But on the other hand, the sum xq + ■ ■ ■ + x n has the distribution 1 + S n . We note that 
1 + S n = 2B n — Ti + l, so that: 

l A l(k)t k = £ l- ] E((x Q + x l + ---+x n ) k )t k = (5.1) 



fc=l,3,... 



fc=l,3, 



M 2Bn _ n+l (t) - M 2Bn -n+i{-t) e* - e~* /e* + e~ l ^ " 



6 Precise formula for the linear part 

Recall (see equation (|4.5|> ) that: 

LIN n (T) = LIN \T ' - l{> ~ r ' / " 



n 



In this section we will derive the precise formula (|6,2p for LIN n (T). 
We have: 



/() + ••• + /« 



3! (^) 3 

The coefficient of the term of this expansion that depends on Ya=o /* linearly is: 



or in terms of the differentiation operator d: 
LIN n {T) 



T'(0) r®(o)_ Z^L^ 



(9 9 3 S 5 



3!(y^) 3 ' 



3!(^) E 



(T)(0). 



Using the formula (|5.ip . we get: 



LIN n {e^ x ) 



•1) 



or in general: 

LIN n (T) 



_ e —d/y/n\ ( e 8/^/n i g — d/y/ri 



(T)(0). (6.2) 
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7 Solution of the extremal problem 

Recall (see remark 14. 4p that we are looking for an odd function T (introduced in 
theorem 14. 7p satisfying: 

sup \T(x)\ = 1, 

— oo<x<oo 

for which the quantity: 

lim sup \fn ■ LIN n (T) 

is maximal. 

The following theorem proves that under some additional conditions on T this ex- 
tremal problem has a solution and that solution is: 

T(x) = sin(x). 

Theorem 7.1. // the function T satisfies the following conditions: 

a) the support of the Fourier transform of T is countable, and the sequence {ujj} c ^ =l is 
linearly independent over the set Z of all integers: 

oo 

T{x) = Y J C 3 e-^ x , (7.2) 
i=i 

b) T is an odd function, 

c) \ \T\\oo = 1 
then 

lim sup • LIN n (T)<—= (7.3) 
and the maximum is obtained for the function: 

T(x) = sin(a;). 

Remark 7.4. In theorem \7.1\ by LIN n (T) we understand the operator LIN n formally 
applied to the series (|7.2p by formula (|6.ip . 



Remark 7.5. In theorem \7.1\ the inequality (|7.3p is understood in the sense that if the 
left hand side is real then the inequality is satisfied. 

Proof. Write 

oo 
3=1 

Note that according to [7] (page 183, paragraph 6.9, theorem 9.3) we have: 

oo 

sup \T(x)\ =^2\cj\. 

—oo<x<oo . . 
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Hence: 

oo 

£M = i. 

3=1 

On the other hand, according to the remark T7.4I we get: 
lim sup y/n ■ LIN n (T) = lim \fn ■ LIN n (T) 



de d2 / 2 



(T)(0). 



We calculate: 

'de*' 2 

Plugging in x = 0, we get: 



(T)(x) 



-/ > ,,:,e-^ /2 e~ ix ^ - Cj . 



3=1 



3=1 

In order to continue with the proof note the following remark: 
Remark 7.6. We have: 

-uj 2 /2 1 

max uoe 1 = —=, 

— oo<w<oo yj e. 

and the maximum is obtained at the point to = 1. Similarly: 



mm ooe 

-oo<o;<oo 



and the minimum is obtained at the point oj = — 1. 
Thus, in order that the maximum of the quantity: 



3=1 



-«?/2 



subject to the restriction: 

oo 

is obtained we need the sequence {uij} to consist of only 2 points, namely: 1 and — 1. 
As a consequence we get that the odd function T that maximizes the quantity: 

'de 92 ' 2 } (T)(0) 

subject to the conditions stated in the theorem is: 



T(x) 



e~ lx + 



and in that case: 



2 

de 92 / 2 



2 

(0(0) 



sin(x), 



1 



□ 
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8 Derivation of the best constant 



Theorem 8.1. We have 



lim inf inf , 

N-+00 VlnN tt(P)=iV Jo 



1 /•! 



\D P \ > 



and 



1 ^ f 1 

lim sup . inf / / 

jv^oo VlnTV tt(^)=A f 7o 7o 



|Dp| > 



256Veln2 
1 

64Veln2 



Proof. Taking T(aj) = sin(x) we get: 



= sm 



n 

J=0 



cos I — ^= I + i ■ fj ■ sin 



cos 



COS 



_n-2 



sm 



£)( E /*/*&)+■- 

V 7 \0<il<j 2 <i3 / 



So that: 



1 f 1 n ■ ffo + "' + fn 
/ Dp sm 

o Jo 



cos — = sm 



) ( ^ 



^0 



£>p/i 



-3! cos 



n-2 



sm 



E 

v 0<ji<j2<73 

For the first term using lemma [2T9l we estimate: 



J! Jo Dpfjjjjj *J + 



cos I — p= I sm 

In 



/n) £ 



JO 



Dpfi > 



1 \ / 1 \ /V2 _2n / 1 \ N2~ 2n 

2 ™" I is ' +0(1) ("+dp"-^>— a b +o(1) ^-"J— • 



J_ 

7^ 



Let us point out that the appearance of: 



here is not surprising as it corresponds to LIN n (sin) (see section [7]) . 
For the second term we apply lemma [2.11I to get: 



3! cos 



n-2 



sm 

n ) V \/n 



fa I 



0<jl<j2<j3 



1 /•! 



JO 



D Pfjlfj2fj3 



< 



14 



1 N2~ n ^ , 1 N2~ n 



< 3I * y 2 n-n < 3, - . y 2 - dd{n _ d) = U-\, 



0<?i<72<.?3<« d=l 

and the higher order terms are estimated similarly. 
Hence we can write: 



Thus since: 
then 

and 

in fact: 



inf HApIIi > inf \D P H n \ > 

2™- 1 <2$(P)=2N<2™ 2"- 1 <2$(P)=2N<2™ 

/ 1 \ N2~ 2n ( 1 

n^ +o(i) )^ (2 "- Ar) — +o b 

2™" 1 < 2Jt(P) = 27V < 2 n 



liminfdjv > 



limsupdAr > 



256y / eln(2) 
1 



64^ In (2) 
1 1 ^ 1 



d 2 n-l > ; : =— 



A /ln(2"" 1 ) 64^/i 6<Veln(2) 



□ 
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